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transition  probabilities  directly,  and  present  a closed-form 
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perform  a gradient  optimization  for  placement  of  stationary 
sea-chers  in  an  ASW  context . 
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PREFACE 

This  report  presents  a new  formulation  in  search  theory  Inspired 
by  Pacific-Sierra  Research  Corporation's  continuing  analysis  of  ASW 
problems,  the  formulation  is  distinguished  by  1)  a realistic  character- 
ization of  the  search  environment  relative  to  prior  analytical  models, 
and  2)  low  computational  costs  compared  to  simulation  techniques.  The 
report  should  be  of  interest  to  operations  researchers  responsible  for 
for  optimising  search  tactics  and  assessing  search  technology. 
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1.0  INTRODUCTION  AND  SUM^\RV 

A fundamental  problem  in  search  theory  is  to  calculate  the  proba- 
bility of  detection  for  searchers  attempting  to  find  a target  whose 
initial  location  and  subsequent  motion  are  characterized  statistically. 
The  searchers'  paths  and  laws  of  detect  ion  are  assumed  known.  The 
solution  of  this  problem  has  important  applications  in  the  optimization 
of  search  tact  ics  and  the  assessment  of  search  technology.  Because  of 
a lack  of  relevant  and  accessible  research  on  the  moving  target  problem, 
operations  analysts  have  often  resorted  to  overly  simple  analytical 
procedures  or  expensive  computer  simulations  to  obtain  numerical  results. 

This  report  presents  u new  formulation  that  provides  an  exact  solu- 
tion foi  the  search  problem,  assuming  a Markovian  search  and  a target 
whose  motion  is  a diffusion  process.  These  assumptions  are  then  relaxed 
to  ai ,'ly  the  search  formulation  to  seveial  important  non-Markovian  tar- 
gets and  searchers.  In  general  our  purpose  Is  to  1)  render  the  search 
formulation  accessible  to  operations  analysts  without  sacrificing  mathe- 
matical rigor;  2)  narrow  the  gap  between  certain  analytical  assumptions 
and  the  behavior  of  actual  targets  and  existing  search  technology;  and 
3)  illustrate  the  utility  of  the  formulation  with  simple  numerical 
examples.  Afcer  a brief  review  of  previous  related  research,  the  pre- 
sent section  summarizes  the  developments  in  this  report. 

1.1  PRIOR  RESEARCH 

In  his  pioneering  work  [1],  Bernard  Koopoan  solves  problems  in- 
volving he  search  for  moving  targets  with  bis  classical  "random  search" 
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detection  model.  Although  it  remains  a viable  analytical  tool  thirty 
years  after  its  introduction,  the  random  search  model  is  limited  by 
omission  of  certain  operational  considerations  such  as  a realistic  repre- 
sentation of  either  target,  or  searche ; motion,  or  an  aeconroodat ion  of 
general  detection  laws. 

Later  results  in  search  theory  have  generally  emphasized  the  optimal 
allocation  of  search  effort  for  detecting  targets  that  are  either  sta- 
tionary or  have  simple  .inds  of  motion.  The  recent  book  by  Stone  |2] 
presents  a unified  treatment  for  much  of  this  wotk. 

Heilman  [3]  provides  a general  formulation  of  the  Markovian  search 
problem  that  is  similar  in  many  respects  to  developments  given  here. 
However,  the  level  of  presentation  and  the  use  of  a Fokker-Planck  equa- 
tion to  describe  target  motion  may  have  limited  the  application  of 
Heilman's  work  by  those  seeking  computational  results  for  actual  search 
problems.  For  comparison,  we  describe  the  work  of  Koopman  and  Heilman 
In  somewhat  more  detail  atter  developing  our  formulation  of  the  starch 
problem. 

1 . 2 SUMMARY 

Section  2.0  considers  the  search  for  a Markovian  target  that  moves 
discretely  in  both  time  and  space.  That  is,  the  target  may  be  situated 
In  any  one  of  a countably  Infinite  number  of  locations,  and  is  capable 
of  jumping  to  another  location  at  each  time  step.  The  initial  distri- 
bution and  the  locat ion-to-locat  ion  transition  probabilities  for  the 
target  are  known,  as  well  as  the  probability  that  the  searchers  will 
find  the  target,  given  its  location.  The  latter  is  determined  by  t he 
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searchers'  pachi  and  lavs  of  letection.  We  aerlve  a recu.s.ve  difference 
equal  ion  for  the  distribution  of  the  target  a*  any  time  conditioned  on 
an  unsuccessful  search  until  that  time.  Thia  . xpression  is  then  used 
to  obtain  the  probability  oi  detection.  Although  it  only  crudely  repre- 
sents the  movement  of  real  targets,  the  discrete  formulation  if.  a good 
model  for  certain  generalized  search  problems,  as  we  show  in  an  example 
Involving  the  Interception  and  decoding  of  clandestine  messages. 

The  discrete  search  problem  of  Sec.  2.0  also  Introduces  the  con- 
tinuous case  treated  in  Sec  3.0.  Working  similarly  but  at  a somewhat 
higher  mathematical  level,  that  section  derives  a nonlinear  integrodif- 
ferential  equation  for  the  target  location  density  conditioned  on  an 
unsuccessful  search.  The  target  i^  assumed  to  move  ae  a diffusion  pro- 
cess with  known  initial  density  and  transition  probabilities.  The 
searchers'  ability  to  find  the  target  is  represented  by  the  Markovian 
search  densirv  introduced  by  Koopman  end  used  throughout  the  literature. 
By  working  with  the  lolnt  events — target  location  and  unsuccessful 
search — we  linearize  rhe  integrodif ferent  ia  1 equation  and  use  its  solu- 
tion to  estahllrh  a simple  expression  for  the  probability  cf  letection. 
Section  3.0  concludes  with  a physical  interpretation  of  the  continuous 
search  formulation,  comparing  it  with  the  treatment  by  Heilman. 

Section  4.0  applies  the  continuous  search  formulation  to  certain 
operational  aspects  of  actual  searches.  The  search  density  of  Sec.  3.0 
is  expressed  as  an  explicit  function  of  the  individual  laws  of  detection 
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ar.d  paths  of  an  arbitrary  number  of  moving  or  stationary  searchers. 
Using  K^'-pcnan's  lateral  range  curve,  we  find  an  approximation  that 
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extends  the  Mr  opt*  ot  t hr  search  formulation  to  many  Impi'i  t ant  no  i- 
Markovian  search  devices.  Thr  approximation  involvrs  the  I’lvoroloti 
o‘  an  Abel  int rural  equal ton  to  obtain  thr  "rqui valrnt " Markovian  law 
of  tit*  tor  l ion  from  attv  wr  1 1 - behaved  lateral  range  curve. 

That  sort  ion  also  extends  the  search  formulation  to  Include  a simnie 
but  important  class  of  non-Markov i an  targets  of  which  the  classical  tier- 
ing datum  is  a special  case.  The  fleeing  datum  Is  a target  that  initially 
Selects  a heading  at  random  in  the  Interval  0°  to  360°,  and  then  move 
on  a straight-line  course  at  a known  constant  speed.  Closed-form  solu- 
tions of  the  differential  equations  derived  in  See.  3.0  are  obtained  for 
such  a target.  These  solutions  are  used  in  a numerical  example  to 
illustrate  the  dependence  of  the  sear*  h path  on  detection  probability; 
the  result  is  computed  to  one  obtained  with  the  pat h- Invar iant  random 
search  model  of  Koopmuu. 

The  last  section  concentrates  on  the  numerical  optimization  of 
search  effort  in  a problem  dealing  with  stationary  sensors  used  to 
detect  a fleeing  datum  with  an  initial  circular  normal  distribution. 

The  sensors  are  assumed  to  have  : definite-range  law  of  detection.  We 
use  a simple  and  Inexpensive  computational  procedure  to  t Ind  the  optimal 
circular  sensor  pattern  that  maximizes  the  probability  of  detection  for 
a search  of  fixed  duration.  The  optimization  is  performed  over  a range 
of  input  paramo  .s  relevant  to  a submarine  search.  For  a specific  set 
of  inputs,  the  best  circular  pattern  with  e gilt  sensors  is  shown  to  he 
superior  to  an  optimal  square  pattern  using  nine  sensors. 
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Finally,  we  consider  ? tie  eit'ect  or  target  course  changes  subsequent 
to  the  initial  heading  selection  of  the  fleeing  (iuUai.  The  search  is 
undertaken  by  sensors  in  a circular  pattern  previously  optimized  fcr  a 
specific  fleeing  datur.  problem.  All  inputs  remain  as  in  that  problem, 
cxcent  that  the  target  chooses  a new  course  every  fixed  time  increment. 
The  surprising  result  is  rhai  additional  course  selections  hardly  in- 
crease the  target's  chance  of  escape;  in  fact,  they  decrease  its  chance 
of  escape  if  not  sufficiently  frequent.  Section  5.0  concludes  with 
suggestions  for  broadening  our  computational  experience  and  extending 
the  analytical  development  of  the  search  formulation. 
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2.0  SEARCH  IN  DISCRETE  SPACE  AND  TIME 

This  section  considers  the  search  for  a Markovian  target  that  moves 
in  discrete  time  in  a space  consisting  of  a denumerable  set  of  states  or 
locations.  The  discrete  search  problem  provides  a structure  in  which 
we  can  consider  the  more  complex  situation  where  the  search  process  and 
target  movement  occur  in  continuous  time  and  space  (i.e.,  where  the  tar- 
get moves  as  a diffusion  process).  The  continuous  case,  treated  in  Sec. 
3.0,  is  the  most  relevant  to  reaL-world  search  problems,  and,  happily, 
will  lead  to  closed-form  solutions  for  the  cumulative  probability  of 
detect  on  in  many  problems  of  interest. 

The  discrete  search  problem,  however,  has  intrinsic  importance  be- 
yond merely  introducing  the  continuous  case  as  demonstrated  by  the  example 
at  the  end  cf  this  section.  In  addition,  search  problems  are  often  part 
of  a larger  model  that  has  a previously  defined  discrete  space  for  target 
motion.  The  discrete  formulation  is  therefore  especially  useful  in  a 
large-scale  simulation  of  a military  operation  where  search  is  one  of 
many  functions  that  contribute  to  its  overall  success.  Such  simulations 
are  handled  almost  exclusively  by  digital  computers  where  friendly  and 
enemy  forces  are  constrained  to  move  on  a finite  grid  in  discrete  time 
steps . 

Ignoring  for  the  moment  the  distinctions  between  the  discrete  and 
continuous  cases,  we  can  make  the  following  informal  description  of  the 
general  cla^s  of  search  problems  addressed:  Find  an  expression  for  the 

cumulative  probability  of  detection  lor  a given  number  of  moving  or 
stationary  searchers  attempting  to  detect  a target  whose  location  t 
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t ■ 0 and  notion  for  t >_  0 are  governed  by  a specified  initial  distribu- 
tlon  and  transition  probability  function,  respectively.  Each  searcher's 
path  (fixed  point  for  a stationary  searcher)  and  law  of  detection  are 
also  given.  Our  approach  Is  tc  first  derive  sn  expression  for  the  location 
density  of  the  target  for  any  tine  t _>  0,  conditioned  on  an  unsuccessful 
search  in  the  Interval  f0,t>.  The  solution  of  this  expression  Is  <-hen 
used  to  establish  P(t),  the  probability  of  detection  in  the  interval  [0,t). 
For  the  discrete  case,  this  procedure  entails  finding  an  expression  for 
P(n),  the  probability  of  detection  by  tine  n,  after  establishing  a re- 
cursion for  the  distribution  of  the  carget  at  n given  that  rhe  target  has 
not  been  detected. 

2.1  laftGET  MOTION 

Let  the  Vocation  of  the  target  at  tine  n be  represented  by 
{*tt,n  «*  0,1,2,...},  an  irreducible  Markov  process  defined  on  the  state 
space  X taken  to  be  the  set  of  positive  Integers  x - 1,2,?,....  Hence, 
target  notion  is  considered  to  be  a discrete  Markov  process  with  respect 
to  both  the  state  variable  x and  the  tine  variable  n. 

At  n • 0,  target  location  Is  characterized  by  the  known  initial 
target  location  density  p (x),  defined  by 


P (x)  - *}  , xeX 

o o 


such  chat 


x-l 


\ 


Ter get  aavtmnt  is  represented  by  the  known  stationery  one-step  tran- 
sition probabi-  li  ties 


*lj  1 VU 

such  that,  for  l.jeX*  and  n * 0.1.2,  — , we  have  as  usual 

(D  1 0 

and 

<u)  Xaj-1 

J-i 

n ordarlng  target  transition  and  search  events,  a prerequisite  for 
the  analysis  to  follow,  we  assuae  that  the  r.L‘‘  target  transition  occurs 
Instantaneously  at  tine  n - t,  n ■ 1,2,...,  where  t <<  ] is  arbitrarily 

snail . 


2.2  SEARCH  PROCESS 

Consider  the  detection  event  defined  by 

B : {target  detection  at  tine  n + e)  , 
n 

where  it  is  assuned  that  all  xt=X  are  searched  instantaneously  at  each 
tine  n + e,  n ■ 0,1, *2,...,  and  (as  in  the  case  of  target  transitions) 
e <<  1 is  arbitrarily  snail.  The  ability  of  the  searcher  to  detect  the 
target  is  characterized  through  the  known  search  density  y(x,n)  defined 
by 


9 


rU,n)  -#{1  | X -x) 
n n 

for  *11  xeX  a®*1  n " 0,1,2,.... 

We  remark  Chat  the  aearch  density  y(s,n)  results  rro»  the  combined 
effect  of  • finite  number  o£  stationary  or  moving  searchers,  each  with 
Its  own  (possibly  unique)  detection  rule  and  search  path  (a  fixed  point 
In  X for  a stationary  searcher).  Furthermore,  the  search — more  generally, 
each  searcher — Key  commence  at  some  time  n'  after  the  Initial  target 
"fix"  with  associated  density  0Q(x);  in  this  case,  y(x,n)  = 0 for  all 
x€£i  and  n ■ 0,1,2, ... ,n'  - 1.  Both  of  thaee  points  are  elaborated 
In  Sec.  4.0,  which  considers  the  applications  for  the  continuous  search 
formulation. 

We  limit  the  present  analysis  to  x Markovian  eearch  in  the  sense 
that  the  result  of  a search  at  time  n + e depenf.a  only  on  the  location 
of  the  target  at  that  tine,  and  not  the  results  of  prior  searches.  More 
precisely, 

*{BJ  V1-  Vr  V2 V 

- y(i,n) 

The  Markovian  property  has  two  interpretations.  If  we  define  the 

search  as  finished  at  the  time  of  the  first  detection,  then  a search  at 

time  n + e implies  the  realization  of  the  events  B , , B B , 

r n-1  n-2  o 

and  the  property  has  a trivial  meaning.  On  the  other  hand,  if  we  choose 
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to  continue  the  search  Indefinitely,  then  the  Markovian  property  allows 
ua  to  count  detections  without  altering  our  search  strategy  (remember- 
ing that  y(x,n)  ia  fixed  a priori).  In  the  development  below,  we  solve 
for  the  probability  of  the  first  detection  by  time  n.  Thus,  both  inter- 
pretations are  viable. 

2.3  INDEPENDENCE  OF  TARGET  MOTION  AND  SEARCH 

On  intuitive  grounds,  the  target  is  unaware  of  the  search,  and  the 
searcher'o  strategy  ia  fixed  a priori;  target  motion  and  search  should 
thus  be  indepen>  ent  in  some  sense.  To  embody  this  intuitive  notion 
concretely,  we  introduce  the  independence  property: 


#{X  .,  •»  1 I X - i. 

utx  n 


B . B B 1 

n n- 1 o 


-<P(X 


n+1 


1} 


This  property  simply  states  that  the  (n+1)  transition  is  independent 
of  the  previous  n searches.  Rewriting  it  in  a slightly  different 
form. 


<*{Xn+l  - J-  Bn 


V1’  Vl’  *n-2’ 


’V 


“^Xn+l“J  I X-“1}^{D* 


X -1} 
n 


■ Y(i,n) 


where  we  have  employed  the  Markovian  search  property  as  well . 
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2.4  SEARCH  THEOREMS 

In  order  to  establish  P(n),  the  probability  of  detection  by  time  n, 
we  introduce  p(x,n),  the  target  location  density  at  time  n,  given  an 
unsuccessful  search  until  that  time — or,  more  simply,  the  conditional 
density.  Thus, 

P(x,n)  «d^(X  -x|  B B } (2.1) 

n n-l  o 


for  xeX«  The  following  theorem  obtains  a recursion  for  p(x,n)  from  the 
givens  pQ(x),  and  y(x,n). 


THEOREM  2.1: 


mL  J.*  ^ * * 3 • a ✓ • \ » 

xr*  Ltvruxi'OUjriui'  Ut'fMJbUy  pU*Uj  HU  UZH  / LttH  LM£  rtzClWHZ>V*i  ejUUX  VOTl 


^U'lfU.n)]  ^ p(i,n) 

P v j ,n+l)  - , n “ 0,  1,  2 (2.2) 


^ ^ [l-Y(l.n)]  p (i,n) 
i-1 


uilh  initial  condition  p(i,o)  - pQ(i),  ieX- 


PROOF: 


From  Eq.  (2.1)  we  have 


P<j,n+1)  -^Un+1‘  i I Bn Bo) 


B_„ "o' 


,/fo' 


n 


B } 

o 


» « * • » 
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^Xn-H  Bn  * Vi Bo} 

^{®n  I Vl ®o} 


UVl •„>  “JVl V 


- IX  1 Xn  ' l-  B»-l B°)^tX"  ' ' 1 B"-‘ 


_ \ &I  n I v mi  R R \ n t < T*  \ 

/ < r.  ' n n-1  c 


From  the  Markovian  starch  property, 


Vl V 


- 1 - v(i.n)  , 


so  that 


i I B B } - / . ll-v(l.n)]  p(l,i 

n n-1  o m 


i1" 
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Similarly, 


^{Xn+l“  ,f  ®n  I Vl* ®o} 


,^{X  - 1,  X - J,  B | B 8 } 

1 n n+1  n 1 n-1  o 


»„  i v *•  1 1 »„-i »„> 

i-i 

m 

•5^‘Vl-J-  *J  V1-  Vl 


from  the  Independence  property. 


Therefore, 


W1-  Vi ‘o' 


•'‘Vl")  V1' 


♦ij  (l-Y(i.n)] 


•'‘Vi')'  »JVi V 


[1-y  (i,n) ] o(i,n) 


which  icnpletes  the  proof. 
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Nov  ve  let  Q(n)  be  the  probability  of  an  unaucceaaful  search  until 
time  n.  Thua, 


«•>  -^Vl'  B„-2 “o1 


(2.3) 


and  the  probability  of  detection  by  n is 


P(n)  - 1 - Q(n) 


(2.4) 


The  following  theorem  establishes  P(n)  from  knowledge  of  p(i,n). 
THEOREM  2.2: 

The  probability  of  detection  by  tine  n ia  given  by 


n-1  — 


- i - 


■ I v 

11  Zj 

n'-O  i-1 


{1-y p(i.n’) 


i •%  c\ 
U*  J9 


where  y(i,o)  - 0 and  p(l,o)  - Pq Cl)  for  all  i€X- 


PROOF: 


From  Eq.  (2.3)  we  have 


\ .^iB  , . • - , B } 
Q(n-fl ) n o 

Q(n)  “ 


d*B  ......  B } 

n-1  o 


"o’ 


\ 

\ 
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- 2>» 


i. 


B 

n 


n-1 


OO 


1-1 


X 

« 


1}.*»{X  - 1 I B B } 

n 1 n- l o 


[l-Y(lfn)J  p (i,n) 


Thus, 


uin) 


n-1  « 

r~i  \ 
i 1 L 

n’-O  i- 


li-y(l,n:)j  n(t,n') 


which  together  with  Eq.  (2.4)  completes  the  proof. 

2.5  DISCRETE  SEARCH  PROBLEM 

As  sn  application  of  the  discrete  search  formulation,  consider  the 
following  generalized  search  problem.  Each  day  Red  selects  a broadcast 
channel  to  send  an  encoded  message  to  remote  operatives.  In  an  attempt 
to  intercept  and  decode  Red's  message,  Blue  also  selects  a channel  each 
day.  Both  Red  and  Blue  are  limited  to  one  selection  dally,  and  we  assume 
that  commensurate  sets  of  channels  are  available  to  each.  We  also 
assume  that  Blue  has  some  knowledge  of  the  process  by  which  Red  aakes 
daily  channel  selections.  The  problem  is  to  assess  Blue's  strategy  for 
Intercepting  and  decoding  Red's  message. 


I {> 


The  "target'1  In  this  case  is  Red's  broadcast  channel.  Let  the 
Markov  process  (X^.n-  1,2,...!  with  state  space  X represent  Red’s  dally 
channel  selection.  Thus,  U X 1*  taken  to  be  the  finite  set  of  integers 
x“l,2,...,l,  then  - x has  the  interpretation  that  Red  chooses  to 
broadcast  on  cltannel  xeX  f°r  day  n.  We  assume  thar,  on  day  zero,  all 

channels  are  equally  likely,  so  p (x)  - I ^ for  all  xGX*  Blue's  knoul- 

o 

edge  of  Red's  broadcast  sequence  is  embodied  by  the  one-step  transition 
probabilities  ^ for  all  i,.]€X- 

Blue ' a strategy,  fixed  a priori , consists  of  the  monitored  sequence 
of  channels  denoted  by  <z^>,  such  that  X fot*  n "0,1,2,....  If  Blue 

has  correctly  selected  Red's  broadcast  channel  on  any  day,  then  he  has 
probability  yq  of  decoding  Red's  message.  Blue's  "search  density"  is 
thus  given  by 

1Y  x - z 

o n 

0 otherwise 

for  all  xeX  and  n*0,l,2t....  The  probabil  l ty  that  by  day  n Blue  will 
have  successfully  decoded  a message  sent  by  Red  is,  from  Eq.  (2.5), 

n-1  I 

P(n)  - 1 - | | ^ *]  1 1 - Y(i,n')  ]p(l,n’) 

n'-O  i-1 


where  p(i,n)  is  found  from  the  recursive  equation  (2.2)  t 
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1 

^ , [ I - rO  .»)  |^{  P ( 1 ,n) 

p(|.n+l)  - l-r  1 - 

k 

^ ^ ( i - y(l,n) ]p(t,n) 

1-1 

with  lnltl.il  condition  p(i,0)  - I ' for  aLl  1GX- 


IH 


0 JlEAR ^H_. JPOR  _A  TARGET  WHOSK  MOTION  IS  A DIFFUSION  PROCKSS 

This  sort  ion  derives  tin*  key  result  of  our  analysts* -a  I inoar  dif- 
ferential equation  from  which  a simple  expression  for  the  probability 
of  detection  is  obtained  when  the  search  is  undertaken  in  continuous 
space  and  time.  The  target  is  assumed  to  move  as  a diffusion  process 
with  a known  transition  probability  function  and  initial  distribution. 

The  search  process  is  the  continuous  space  and  t imr  analogue  of  the 
discrete  Markovian  search  discussed  in  Sec.  2.0. 

The  importance  of  the  continuous  search  problem  is  twofold.  First, 
physical  considerations  suggest  a continuous  space  and  time  representa- 
tion of  target  and  searthor  behavior.  Real  targers  generally  cannot 
move  over  finite  distances  in  zero  time,  and  existing  search  devices 
generally  operate  eont  inuously  in  tl.ae.  Happily,  these  observations 
suggest  that  the  processes  used  to  define  the  search  problem  have  little 
difficulty  meeting  the  regularity  and  continuity  conditions  required  for 
the  derivations  below. 

Second,  the  continuous  search  formulation  involves  expressions  that 
are  easier  to  manipulate  than  the  recursion  of  Theorem  2.1  for  tin*  dis- 
crete case.  Indeed,  Sec.  4.0  presents  a closed-form  soluti*  n for  a 
classical  real-world  problem,  the  search  for  a fleeing  datum.  Only 
sltnpLe  numerical  integrations  are  required  to  compute  he  (reliability  of 
detection  for  this  problem. 

The  analysis  of  this  section  is  r.ecessurtly  carried  out  at  a signi- 
ficantly higher  mat Semat lea  1 level  than  that  of  Sec.  2.0  for  the  discrete 
search  problem.  However,  the  structure*  of  the  derivations  remains  the 
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savxie:  establish  a rela<  ion  tor  the  conditional  target  location  density, 

and  derive  an  expression  tor  the  probability  of  detection  as  a function  of 
the  co  nditional  density.  Th«  li»i:egrodi  fferential  equation  for  the  con- 
ditional density  is  linearized  by  considering  the  joint  density  (i  e., 
the  density  associated  with  the  Joint  event,  target  location,  and  un- 
successful search).  Section  4.C  solves  the  linear  integrodif ferential 
equation  for  the  joint,  density  under  the  assumption  of  the  classical 
flee Ing-da turn  target  motion  and  an  arbitrary  Markovian  search. 


3.1  TARGET  MOTION 

Consider  the  Karkov  process  {X(t),  t >0 ) with  continuous  time  p ra- 
meter  t«=  (0,®)  and  state  space  X taken  to  be  E^,  the  two-dimension*. 
Euclidean  space.  The  location  of  the  target  at  time  t > 0 is  represented 
by  the  random  variable  X(t)eXwith  coordinates  (X^Ct),  X^t)),  - <®  < 

-<i  (t)  < <*. 

At  ' w 0,  the  target  is  located  according  to  the  known  iniriaJ. 
density  p (x) . Thus, 

Cj 


and 


p (x)dx  K^{X(0)G  dxl 
o 


(3.1) 


where  dx  is  an  infinitesimal  element  of  area  containing  the  point  xGX  • 

2 

and  the  integral  is  over  E (as  will  be  the  caco  for  all  integrals  dcIow, 
unless  otherwise  noted). 

The  known  transition  probabilities  of  (X(t),  t>G}  are  given  by  the 
function  »ji(x,t;y,T),  i > t,  x , yy_X , such  that 


20 


*(x,t;y,T)dy  - ^(X(t)  €dy  | X(t)  - x)  (3.2) 

and 

j <i(x.t;y,x)dy  - 1 . (3.3) 

Since  (X(t),  t>0}  is  a Markov  process,  it  satisfies  the  Chapman-Kolaogorov 
equation;  namely,  for  se(t,f), 

\Kx,t;y,t)  - J i>(x,  t; z, s)  tKz,s;y,t)dz  . (3.4) 

We  assume  that  ^(x,t;y,t)  has  a derivative  with  respect  to  t at  i ■ t, 
so  that  we  may  write  for  small  At 


£(x,t;y,r.+At) 


S(y-x)  + At 


3t 


+ o(At) 

T-t 


(3.5) 


where  £{•)  denotes  the  two-dimensional  Dirac  delta  function. 

We  will  also  need  a set  of  regularity  conditions  that  limit  the  be- 
havior of  the  target  over  small  intervals  of  time.  These  conditions  will 
be  recognized  as  those  used  in  the  derivation  of  the  Kolmogorov  diffusion 
equations.  Wi.th  taken  to  be  a circle  of  radius  5 > 0 and  center  x, 
the  regularity  conditions  for  x m (x.,x_),  x.yGX.  te{0,«),  and  small  At 
are  given  by 

j «,(x,t;y,t+At)dy  - o(At)  , (3.6) 


J (yi~xi)  iKx, t;y,t+At)dy  - a^x.OAt  + o(At)  , 


i - 1,2  , (3.7) 
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S« 

* b^j(x,t)4t  + o('lt)  , i,J*l,2  * (3.8) 

and 

/ (yi~Kl)U  ^yj~XJ>V  ,*'(x*f-5y»t+At)dy  “ °(At)  , l.J-1,2  , (3.9) 

S6 

where,  for  Eq.  (3.9),  u,v  0 and  u + v > 2. 

3.2  SEARCH  PROCESS 

The  search  is  undertaken  continuously  ir  tine  and  is  Mai covian  in 
uatuie.  If  we  denote  the  aececcion  event  for  t > t by 

B(t,x) : (target  detection  within  [t,t)}  , 

then  the  known  search  density  Y(x,t)  is  defined  by 

y(x,  t)At  + o (At)  - ;j»{B(t,t+At)  | X-x)  (3.10) 

for  x€X.  te[0,«),  and  small  At.  Note  that  the  probability  in  Eq. 

(3.10)  is  conditioned  on  a stationary  target  at  x.  We  will  expend  con- 
siderable effort  to  suit  our  definition  of  y(x,t)  to  the  case  of  a 
target  that,  while  moving,  satisfies  'he  regularity  conditions  (3.6)  to 
(3.9).  To  this  end,  we  require  that,  for  all  xeX*  te[0,«>),  a constant 
M exist  such  that 

y(x,t)  <_  M < «®  ; (3.11) 


and,  for  all  x,y<=X»  te(0,»),  and  small  At, 
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1 X(t«t)  - y,  X(t)-x> 


- |y(x,t)  + V 1/n!  Tdj^a/ax.  + 1 Ry (*, 

( n-1,2  J X"y 


t) 


+ o(|y-xj  ) 


At  + o(At) 


where  we  have  made  use  of  the  notation 


(1)  ^ - (.Vj-Xj  ) 

(it)  [d13/3x1  + d23/3x2]2 

- [d232/3Xj  + 2dld?a2/3x13x2  + d232/3x2] 
(iii)  o ( i y-x  j £ ) - o(d^)  + o(d],d2)  + o(d2) 


(3.12) 


Furthermore,  we  assume  that  the  first-  and  second-order  spatial  partials 
of  y(x,t)  exist  and  are  bounded  for  ail  xGX  and  t£[0,®). 

Equation  (3.12)  bears  further  explanation.  It  embodies  the  notion 
th  t,  if  we  know  that  the  target  remains  in  the  immediate  neighborhood 
of  x over  a small  interval  of  time  At,  then  we  can  express  the  proba- 
bility of  detection  in  [t,t+At)  by  an  expansion  of  y(*,t)At  + o(At) 
about  the  point  x.  Indeed,  Lemma  3.1  below  proves  that,  since  the  target 
cannot  leave  the  neighborhood  of  x in  a small  interval  of  time,  it  is  un- 
necessary to  condition  on  the  nearby  endpoint  X(t+At)  * y. 

As  in  the  discrete  case,  we  muse  formalize  the  Markovian  nature  of 
the  search.  We  thus  require  that  the  search  process  satisfy 
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^{B(t,T)  | X(t)  * x,  ¥(0,t)}  - | X(t)  ■ x)  (3.13) 

and 

i?{B(t,t+At)  | X(t+At)  ■ y,  B(0,t)}  - ,^{B(t,t+At)  |X(c+At)-y}  + o(At) 

(3.14) 

for  all  x,y<£X»  te  10»*),  and  small  At.  Finally,  the  independence  of 
search  and  target  motion  is  embodied  by  the  obvious  requirement  that 

,^{X(r)€dy  | X(t)  -x,  B(0,t)}  - ^{X(t)6dy  | X(t)  « x)  (3.15) 

for  all  x.yeXi  t > t,  and  t,te|0,»). 

3. 3 CONDITIONAL  DENSITY 

The  conditional  density  p(x,t)  is  defined  for  dx,  an  inf lnltesimal 
element  of  area  containing  x,  by 

P(x,t)dx  - ^{X(t)6dx  j B(0,t)>  , (3.16) 

where  xGX  and  tG(0,»).  This  subsection  derives  a nonlinear  integro- 
dif ferential  equation  for  p(x,t)  in  terms  of  the  givens  PQ(*)*  «Kx,t;y,T), 
and  y(x,t). 

The  following  lemma  extends  the  definition  of  y(x,t)  to  a moving 
target  that  satisfies  the  regularity  conditions  given  by  Eqs.  (3.6)  to 
(3.9). 


LEMMA  3.1: 

Given  the  process  (X(t),  t>0)  representing  target  location , 
the  search  density  y(x, t'  satisfies  the  following  conditions  for  snail  At.: 

\ 

\ 

\ 
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^(B(t,t+At)  | X(t)  - x}  - 1 - y(x,t)At  4-  o(At)  (3.17) 

and 

B.  i?{8(t,t+At)  | X(t+At)  - y}  - 1 - y(y,t)At  + o(At)  . (3.18) 


PROOF: 

For  par*;  A,  we  write 
^(B(t,t+At)  j X(t)  «x} 

“ /^(B(t,t+At)  | X(t+At)  **  y,  X(t)-x)  d*(X(t+At)  e dy  | X(t)  - x} 
which  by  Eqs.  (3.2),  (3.3),  and  (3.12)  becomea 


»ftj  f s.  a.i,Aa.\  I * 

I Avty  - 


- 1 

^ j 


- 1 - y(x,t)At 


-e-  d23/3x2j  Y(x,t)  + o( |y-x I <')| 


x-y 


4»(x,t;y,t+At)dy  + o(At) 


The  integral  in  this  last  expression  need  be  performed  only  over  S^,  since 
the  iirat-  and  second-order  partials  of  y(x,t)  are  bounded  and  Eq.  (3.6) 
applies.  Furthermore,  Eqs.  (3.7)  to  (3.9)  imply  that  performing  the  inte- 
gral over  Sg  will  result  only  in  At  and  o(At)  terms  proving  part  A of  the 
lemma.  The  proof  of  part  B is  similar,  if  we  consider  the  process 
(X(t) , t>0}  reversed  in  time  and  impose  similar  regularity  conditions. 

We  will  need  the  following  lemmas  tc  obtain  our  main  result  in 


Theorem  3.1. 
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LEMMA  3.2: 


i?{B(t,t-»*t)  | B(0, t)  } - 1-At  J y(x.t)  p(x,t)dx  + o(At)  . (3.19) 


PROOF: 


Using  Ec9.  (3.13),  (3.16),  and  (3.17),  we  quickly  obtain 


^(B(t,t+Ac)  | 1(0, t)> 

- /^(B(t,t+At)  |x(t)-x)  ^{X(t)cdx  | B(0,t)} 
• J [1  - y(x» t)At  + o(At) ] p(x,t)dx  , 


which  proves  the  leasa. 


LEMMA  3.3: 


For  dy,  ok  infinitesimal  element  of  area  containing  ytX  and 


email  At, 


^»{X(t+At)edy  | B(0,t) } 


- dy  < p(y,t)  + At 


p(x,t)dx  + o(At)|  . (3.20) 


PROOF: 


With  help  from  Eqa.  (3.2),  (3.5),  (3.15).  and  (3.16),  we  obtain 


#{X(t+At)€dy  j B(0,  t)  } 


j^{X(t+At)€dy  | X(t)  - x}  ,/{X(t)£dx  | B(0,t)} 


- dy  j t;y,  t+At.)  p(x,t)dx 


THEOREM  3.1: 

Given  the  inputs  po(x),  i^(x,t;y,T),  and  y(x»t)  for  the  continuous 
search  problem , the  conditional  target  location  density  p(y,t),  yeX> 
te [0,™)  a the  oolutxon  of  the  nonlinear  integrodif ferential  equation 

p (x, t)dx  + p (y, t)  [T (t)  - y(y,t) ] (3.23) 

T«t 

with  initial  condition  p(y,0)  ■ po(y)  for  a ^ yeX- 


I 

I 
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PROOF: 

From  Eq.  (3.16),  wp  have 

p (y, t+At)dy  - ,^{X(t+Ai) G dy  | B(t,c+At),  B(Q.t)} 

_ ,y{X(t+At)G  dy,  B(t,t+At)  |B(0,t)} 

{B(t,C+At)  | B(0, c) ) 

_ ,/{X(t+At)edy  1 B(0,t)}  ;/{B(t,t+At)  | X(t+At)  ■ y.  B(O.t)) 
i/{B(t,t+At)  | B(0,t)} 


wot ing  that,  for  snail  e. 


(1-t)-1  - 1 + t + 0(e) 


wxth  the  help  of  Eq.  (3.72)  we  can  write 


j^{S(t,C+At)  | B(0,t)j  1 


1 + r(t)At  + o (At) 


Together  with  Eq.  (3. 14)  and  Lemmas  3.1  and  3.3,  this  expression  results  in 


P (y , t+At) 


p(y,t)  + At  /[3^V'37',V,~]  P(x,t)dx 

+ o(At)j  [1-  y(y,t)At  + o(At)]  [1  + I (t)  + o(At)] 


p(y.O 


o (x, t)dx 


+ p(y,t)  J T ( t ) - y (y , t) j 


+ o (At) 


Transposing,  dividing  by  At.  and  letting  At  ->  0 completes  the  proof. 


i 

< 

1 

t 

i 

i i 


j 


j 
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Corollary  3.1  shows  that  the  solution  of  Eq.  (3.23)  Is  indeed  s 
proper  density  for  all  t€[0,»). 


COROLLARY  3.1: 

The  solution  of  the  equation  for  the  conditional  density  pre- 
sented in  Theorem  3.1  satisfies  J p(y,t)dy  - 1 for  all  t€[0,»). 


PROOF: 


Let  w(t)  * J p(y  t)dy  for  all  te[0,«).  Then  Eq.  (3.23) 


becomes 


3/3t  w(t)  • u(t)  + v(t)  , 


u(t)  " / P (x.t)dxdy 

y x T“c 

- f p(x,t)  [3/3t  f*(x,t;y,i)dy]  dx 
* 1 v JT-t 


r(t)  -Jp(y.t)  (r(t)  - y(y.t)  ]dy 


r(t)  [w(t)  - 1] 


By  Eq.  (3.3),  however,  u(t)  - 0 for  te[Of°°),  so  that 


3/3t  w(t)  - T(t)  [w(t ) - 1 ] 


The  solution  of  this  equation  is 
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t 

/ node 

w(t)  • 1 + K e° 

But  w(o)  * 1 by  Eq.  (3.1),  completing  the  proof. 

We  are  now  prepai  ed  to  consider  the  probability  of  detection  P(t). 
Specifically, 

Q(t)  - ^*{B(0,t)}  (3.24) 

and 

P(t)  - 1 - Q(t)  (3.25) 

for  all  telO,*).  The  following  theorem  obtains  P(t)  from  the  search  in- 
tensity T(t). 

THEOREM  3.2: 

The  probability  of  finding  a target  by  time  t when  the  target 
moves  as  a diffusion  proceed  is  giver. : by 

t 

P(t)  - 1 - exp [-J  T(OdU  , (3.26) 

0 

where  the  search  intensity,  r(t)  « Jy(x, t)  p(x,t)dx,  is  expressed  in  terms 
of  p(x,t),  the  solution  to  the  equation  of  Theorem  3.1. 

PROOF: 

From  Eqs.  (3.22)  and  (3.24), 

Q(t+At)  - ,/{B(0,t-t-At)} 

- ,^{B(0,t)}  ,y{B(t,t+At)  | B(0,O) 


\ 


Q(t)  [1  - r(t)At  + o (At) ] 
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Thus. 


Q(t)At  llCi  At 


and,  after  letting  At  0t  and  noting  that  Q(0)  “ 0,  we  obtain 


L. 

log  Q(t)  - - J r(r)d^ 


proving  the  theorem. 

3.4  LINEAR  SEARCH  EQUATION 

The  nonlinearity  of  Eq.  (3.23)  makes  closed- form  solutions  for  the 
conditional  density  p(x,t)  and  probability  of  detection  P(t)  difficult 

f"  r*  ^V.  Rv  Pnrtft<Hpr<ni>  t-  Km  iWwf  #-  nenn  e 1 4<%nn4  ^«r  n 

~ •“  — j • “/  ~ w “••{)  w««w  v koi  w ax/cuwxuii  vtuuox  k j yut  , um/t  r. 

simply,  the  Joint  density)  p(x,t),  we  can  obtain  a linear  search  equa- 
tion and  a simpler  expression  for  P(t).  Thus,  we  define  for  all  xGX 
and  t e [0,») 


p(x,t)dx  - d^{X(t)t=dx,  B(0,  t)  ) , 


(3.27) 


and  quickly  note  that 


p(x,t)  - p(x,t)  [1  - P(t)] 


(3.28) 


and 


P( 


t)  - 1 - j p(x, t)dx 


(3.29) 


u 


THEOREM  3.3: 

The  joint  target  location  density  p(y,t),  y^X  > teIO,«)  satisfice 
the  linear  integrodifferential  equation 


3P.& 


H 


^(x,t;y,T) 


3- 


p(x, t)dx  - y(y,t)  p(y.t)  , (3.30) 


T-t 


vith  initial  condition  p(y,0)  « P0(y)  for  all  yeX- 


PROOF: 


The  proof  is  to  simply  transform  Eq.  (3.23)  by  Eq.  (3.28). 


Although  the  Joint  density  has  less  intuitive  appeal  than  the  con- 
ditional density,  it  clearly  represents  the  superior  analytical  tool  for 


solving  search  problems.  Equations  (3.29)  and  (3.30)  will  thus  be  ri»e 


operative  equations  as  we  apply  the  results  of  this  section  to  the 
numerical  problems  of  Secs.  4.0  and  5.0. 


3.5  PHYSICAL  INTERPRETATION  OF  THE  SEARCH  FO'MULATKN 


The  nonlinear  equation  for  the  conditional  density,  Eq.  (3.23),  and 
the  linear  equation  for  the  joint  density,  Eq.  (3.30),  have  simple  but 
important  physical  interpretations.  Both  equations  consist  of  two  terms — 
the  diffusion  term  resulting  in  density  changes  due  only  to  target  move- 
ment, and  the  search  term  causing  density  changes  due  only  to  the  effect 
of  search.  For  example,  consider  the  search  for  a stationary  target, 
i.e.,  a target  with  the  trivial  transition  probability  function  4>(x,t;y,i) 
«(y-x  ).  In  this  case,  Eq.  (3.23)  becomes 


^3*^  - - P (x, t)  [ T (t)  - Y (x,  t)  ] 


(3.31) 


\2 


and  p(x,t)  Is  affected  only  by  the  search  density  y(x*t-).  Alternatively, 
with  no  search — i.e..  If  y(x,t)  - 0 for  all  xe  \ and  t*  |<3,"’) — we  obtain 
the  search-free  equation 


t 


3 T 


T) 


0 (x , t ) dx 
t 


(3.32) 


where  p(x,t)  is  affected  only  by  the  motion  of  the  target.  Similar  observa 
tlons  can  be  made  for  the  joint  density,  p(x,r;- 

The  situation  is  illustrated  In  Fig.  3.1,  where  It  is  assumed  that 
the  initial  density  p (x)  was  circular  normal,  and  that  at  t *»  0 the 

O 

target  chose  a heading  from  a uniform  distribution  on  [0,2x)  a d fled  at 
a known  speed  (the  classical  fleeing  datum).  The  (x^,x2)  plane  is  the 
search  region  and  the  x^  axis  plots  p(x,t)  for  t - T > > 0,  where  the 

search  by  a single  .novlng  searcher  commenced  at  . Tnus,  p(x,t)  Is  tne 
solution  to  Eq.  (3.32)  for  t<=[0,T^),  and  Eq.  (3.23)  for  tf"(Tj,T).  The 
details  of  the  solution  are  taken  up  in  Sec.  A.0;  but  for  now  we  note  that 
the  wide  depression  in  the  center  of  the  p(x,T)  .urface  is  caused  by  the 
Bessel  function  spreading  of  a fleeing  datum  density,  as  first  described 
by  Koopman  [1],  and  the  narrow  depression  along  the  ridge  is  due  to  an 
unsuccessful  search  following  this  path. 

Intuitively,  we  expect  that  if  we  search  an  area  unsuccessfully, 
then  the  object  of  our  search  is  less  likely  to  be  located  in  tliat  area. 
This  phenomenon  caused  the  depression  along  the  ridge  of  p(x,T)  in  llg.  3.1 
Analytically,  Eq.  (3.23)  tells  us  that,  for  points  under  the  influence 
of  the  searcher  (i.e.,  those  for  which  r(x,t)  dominates  its  average 
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r(t)  - j y(x, t)  p(x,t)dx),  the  change  of  p(x,t)  with  time  due  to  the  search 
is  negative.  The  effect  of  T(t)  in  Eq.  (3.23)  is  tc  elevate  p(x,t)  for 
points  relatively  free  of  the  influence  of  the  searcher,  resulting  in 
yp(x, t)dx  ' 1 for  all  te[0,“>)  as  shown  in  Corollary  3.1.  (In  the  linear 
Eq.  (3.30),  the  search  term  f - y(x,t)  p(x,t)]  is  always  nonpositive  so 
rhat  j p(x,  t)dx  » 1 - P(t)  <.  1 for  all  te[0,“).) 

Tlie  search-free  equation,  Eq . (3.32),  is  a direct  result  of  the 
Chapman-Kolmogorov  equation,  Eq.  3.4).  Indeed,  we  can  write  Eq.  (3.4) 
for  t £ [0,t ) as 

vO-,0, y,-t)  »/^(  z,t;v.v)  iMx.0:z,  t)dz 


Integrating  over  the  initial  density  p (x)  results  in 


yii-(x,L);y,T)  PQ(r)dx  * f';(  z,  t;y,t)  |^(x,n,zyt)  Pc(x)dxj 


dz 


or,  in  the  absence  of  search, 


p(y,  *')  ™ f If'(z,  t;y,x)  p ( ;: . r ) d i 


so  tha  t 


3t 


, f[a»J £«.» t;y,T)] 

T»t  J l 


o (z,  :)dz 


T-t 


which  is  the  search-free  equation  (3.32). 

Finally,  we  remark  that  the  only  difference  between  Eq.  (3.23)  and 
the  equation  presented  by  Heilman  [ 3J  is  in  the  diffusirn  term,  because 
Heilman  starts  with  a Fokker-Planck  equation  for  <Kx. t;y,r ) , rather  than 
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t|/(x,t;y,t)  Itself.  (Thus,  by  a manipulation  similar  to  the  one  above 
for  the  Chapman-Rolmogorov  equation,  one  could  obtain  a Heilman  search- 
free  equation  from  the  Fokker-Planck  equation.)  The  specification  of  the 
motion  of  an  actual  target  by  those  responsible  for  the  search  would, 
needless  to  say,  iPad  directly  to  a representative  transition  probability 
function,  rather  than  a Fokker-Planck  equation. 


This  section  applies  the  results  of  Sec.  3.0  to  certain  operational 
aspects  of  actual  searches.  The  emphasis  is  on  target  motions  and  search 
procerses  that  ostensibly  violate  the  Markovian  assumptions  needed  to 
derive  Eqs.  (3.23)  and  (3.30),  on  which  th  continuous  search  formulation 
is  based.  This  section  should  be  regarded  as  an  introduction  to  the  work 
necessary  to  bridge  the  gap  between  the  analytical  assumptions  made  for 
tractability,  and  the  behavior  of  actual  targets  and  existing  search 
devices. 

Subsection  4.1  details  the  straightforward  generalization  of  the 

f 1 « 1 ... I ..  — J — _ C J 
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or  moving  searchers.  Recognizing  that  search  devices  are  often  not 
Markovian  in  the  sense  of  F.q.  (3.13),  we  present  a good  "equivalent" 
Markovian  search  density  in  subsection  4.2  by  inverting  the  integral 
equation  for  Koopman's  lateral  range  curve.  Subsection  4.3  extends  the 
search  formulation  to  handle  the  simplest  kind  of  non-Markovian  target 
motion  for  which  the  classical  fleeing  datum  is  a special  case.  Closed- 
form  solutions  of  the  search  equations  are  obtained  tor  a fleeing  datum 
search.  Finally,  to  illustrate  the  generality  of  the  search  formulation, 
subsection  4.4  presents  a numerical  exaap'  in  which  the  probability  of 
detection  is  calculated  as  a function  of  time  for  two  paths  followed  by 
a moving  searcher  attempting  to  find  a fleeing  datum;  this  result  is 
compared  to  one  obtained  with  a descendent  of  the  path-invariant  "random 
search"  model  first  proposed  by  Koopman. 
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4.1  MULTIPLE  SEARCHERS 

The  search  density  y(x,t)  defined  by  Eq.  (3.10)  for  a stationary 
target  is  actually  the  aggregate  effect  of  the  total  search  effort. 

The  precise  manner  in  which  y(x,t)  depends  on  the  activities  of  many 
individual  searchers  does  not  affect  the  developments  of  Sec.  3.0,  but 
is  central  to  the  application  of  the  continuous  search  formulation. 
Thus,  below  we  construct  y(x,t)  from  the  paths  and  laws  of  detection 
for  each  of  an  arbitrary  number  of  moving  or  stationary  searchers. 

Consider  a search  undertaken  by  n searchers  such  that  the  it*1 
searcher  follows  the  known  path  z^tJeX  for  t G [T^T^)  ♦ where 
T^  > 0 is  the  time  at  which  the  iC^  searcher  commences  his  effort,  and 
Ii2  ii.00)  is  the  time  at  which  he  abandons  the  search.  Obviously, 
for  a stationary  searcher,  Zj(t)  • for  t € [T^.T^)  • The  law  of 

detection  for  the  i searcher  is  defined  in  the  same  manner  as  y(x,t); 
i.e.,  the  probability  that  the  i^  searcher  will  detect  in  the  small 
interval  [t,t+At),  given  that  the  target  is  at  xGX»  is 

yi[x,z1(t) ] A t + o (At) 


The  n searchers  are  assumed  independent  in  the  sense  that,  from  Eq.  (3.10) 


U 

nr 

y(x, t)At  + o(At)  »1-|  I {]  - Y1Ex,zi(t) ]At  + o 


-i 

1 

-l 

n 

(At)} 


“ At  <V"S Y.  (x,z,  (t)  ] + o (At) 
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The  search  density  y(x,t)  is  thus  expressed  in  terns  of  the  in- 
dividual laws  of  detection  Yj_l*»*j(t)  ],  i*l,-,...n,  by 


Y(x,t) 


i-l 


Y1lx,xi(t)] 


(4-1) 


Clearly,  in  order  that  the  search  process  continue  to  satisfy  Eqs.  (3.11) 
through  (3.15),  the  laws  of  detection  for  the  individual  searchers  must 
satisfy  similar  conditions. 

Noting  that,  for  all  xeX>  y^fx.z^t)]  “ 0 when  t < or  t > 
we  can  define  Tj  - aa  the  aeccrch  duratron  where 


and 


X1 


l<i<n 


lil 


l<i<n 


12 


Obviously,  T2  > and  y(xvt)  " 0 *or  t < or  t >_  is  commonly 

referred  to  as  the  "time  late"  in  military  applications. 

4.2  NON-MARKOV IAN  SEARCH:  INVERSION  OF  A LATERAL  RANGE  CURVE 

A sensor  commonly  used  to  search  for  submarines  operates  approxi- 
mately as  follows  (4]:  A signal  that  may  indicate  the  presence  of  a sub- 
marine is  integrated  over  a "sliding"  time  interval  of  fixed  length. 

If  the  value  cf  the  integral  ever  exceeds  a predetermined  threshold, 

Chen  the  sensor  is  activated.  Such  a sensor  is  not  Markovian  in  the 
sense  of  subsection  3.2  because  of  the  "memory"  property  of  the  integral 
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operator.  The  present  subsection  develops  a useful  approximation  for 
non-Markovian  search  devices  that  allows  us  to  apply  the  search  formula- 
tion to  many  situations  where  such  devices  are  employed. 

Whether  or  not  a search  device  is  Markovian,  it  can  be  chare  terized 
by  the  so-called  lateral  range  nii'Ve,  q(0*  The  function  q(£)  is  defined 
as  the  probability  of  detecting  a stationary  target  when  the  searcher 
follows  a straight-line  path  of  infinite  length  and  closest  approach  C 
to  the  target  [1].  The  lateral  range  curve  can  be  approximated  by 
empirical  data  or  computed  directly  from  physical  arguments. 

Consider  a search  employing  a single  Markovian  sensor  with  law  of 
detection  Y^[x,z^(t)],  which  is  a function  of  only  the  target-to- 

searcher  Euclidean  distance  r * |x-z^(t)|.  If  the  target  is  stationary 

? 2 
at  the  origin  ot  E , and  the  searcher  is  stationary  at  z ^ * (£,p)6E  , 

then  denote  the  search  density  by  y(r)  ” 7,10,2,].  Now  let  the  searcher 

move  at  speed  w along  the  path  from  (£,  - “>)  to  (£,<=),  thus  obtaining  the 

probability  of  detection  q(£)  from  the  well-known  expression  11] 

q(£)  - 1 - exp  'Ol  , (4.2) 

where 


— 00 


(Note  that  Eqs.  (4.2)  and  (4.3)  can  also  be  obtained  by  use  of  Eqs. 
(3.29)  and  (3.30)  for  iKx.tjy.t)  = 6(x)  and  p(y,o)  * 6(y).)  Clearly, 
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q(£)  is  the  lateral  range  curve  for  a Markovian  device  with  law  of 
detection  fx. z^(t) ] . 

The  following  question  is  reasonable  to  ask  at  this  point:  If  we 

know  the  lateral  range  curve  for  a radially  symmetric  non-Mar kovian 
device,  can  we  calculate  the  "equivalent”  Markovian  law  of  detection  by 
inversion  of  Eq.  (4.3)?  As  shown  below,  we  can  generally  invert  the 
lateral  ra  ge  curve  and  then  use  the  resulting  Markovian  law  of  detection 
as  input  to  the  search  formulation. 

In  order  to  invert  Eq.  (4.3),  write 

F(0  - ^ /,(V?  + ^)dl1 

n 

OO 

ml  f Y(r)rdr 

-y  (t2.ev/2  • 


2 2 1/2  2 2 
where,  as  before,  r ■ (£  +w  ) . Now  let  r ■ 1/s  and  £ ■ 1/a,  so 

that 


1/2  f , -1/2. 

8 (ti-8) 


)ds 

1/2 


or 


where 


G(a) 


A(s)ds 


(a-s) 


1/2 


(4.4) 


G(a) 


w a 


-1/2 


F<a-1/2) 


and 


w . -3/2  , -i/2, 

X(8y  - S Y (*  ) 


Equation  (4.4)  is  an  Abel  integral  equation  [3],  the  solution  of 


which,  for  continuously  differentiable  C(a),  is 


x<.)  - -m- + 1 c vm* 


,in  "J. 


(4.5) 


(s-a) 


As  an  example  of  the  use  of  Eq.  (4.5),  consider  the  lateral  range 
curve  for  the  inverse  cube  law  c f detection  [1] 


q(£)  - 1 - exp 


I wf 

L ' 


where  w is  searcher  speed  and  k is  a constant  dependent  on  several  aspects 
of  the  search  that  need  not  concern  us  now.  Thus, 


F(0  - -^2 
w£ 


and,  from  Eq . (4.5), 


G(a)  • 2k  a 


X(s)  « — 


f da 

a1/2 (s-a) 


S 

_ 2k  {*  du 

n (s-u2) 
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Therefore,  as  expected  y(r)  ■ k/r  , the  search  density  for  the  inverse 
cube  law  of  detection. 

Although  we  have  shewn  how  to  obtain  the  equivalent  law  of  detection 
for  a non-Markovian  search  device,  we  have  not  shown  how  good  an  approxima- 
tion it  is  for  calculating  the  probability  of  detection  using  the  search 
formulation  cf  Sec.  3.0.  Obviously,  it  is  an  exact  approximation  for 
infinitely  long  straight-line  search  paths,  wherein  lies  the  key  to  a 
general  evaluation  of  the  inversion  technique.  If  the  motion  of  the 
searcher  relative  to  the  target  is  nearly  linear  at  almost  constant  speed 
over  the  detection  range  of  the  aevice,  then  the  Inversion  technique 
represents  a good  approximation.  Further  analysis  is  necessary  to  estab- 
lish the  accuracy  of  the  technique  under  more  general  search  conditions. 

4.3  CONDITIONALLY  MARKOVIAN  TARGETS:  THE  FLEEING  DATUM 

The  search  formulation  is  applied  here  to  a member  of  the  simple 
class  of  non-Markovian  target  motion  described  as  follows:  At  t * 0,  a 

realization  of  an  n-dimensional  random  variable  a is  obtained  from  the 
known  density  f(a).  Target  motion  for  t > 0 is  Markovian  with  con- 
ditional transition  probability  function  ^(x, t;y,r) . Let  the  solution 

of  Eq.  (3.30)  for  (x,t;y,r)  be  denoted  by  p (x,t).  The  joint  density 

a ot 

r 

is  then  given  by  the  n-dimensional  integral  p(x,t)  * J p^(x, t ) f (a)da , 
and  the  probability  of  detection  remains  P(t)  - 1 -f p(x,t)dx.  This 


class  of  motion,  which  we  shall  call  conditionally  Markovian,  is  a 
generalization  of  Stone's  conditionally  deterministic  motion  [2),  where 
knowledge  of  a at  t » 0 implies  deterministic  motion  for  t 5 0.  The 
fleeing  datum,  which  we  now  treat  in  detail,  is  an  important  example  of 
conditionally  deterministic  target  motion. 

The  fleeing  datum  moves  at  a known  constant  speed  u on  a straight- 
line  path  after  chooding  a heading  0 at  t * 0 from  a uniform  distribu- 
tion on  £ 0 , 2 n ) . If  0 is  measured  counterclockwise  with  respect  to  the 
coordinate,  then  the  conditional  transition  probability  function  for 
the  fleeing  datum  is  given  by  the  two-dimensional  Dirac  delta  function 


,!,„fx,t:Y  T • “ -vl  — \r  it  - n i 

T Q '■  > *»  / i */  K / ’Ov‘  •'/  i 


where  » u(cos  0,  sin  0)  is  the  random  velocity  vector  with  magnitude 
u and  direction  ft.  Substituting  Eq.  (4.b)  into  Eq.  (T. 30)  results  in 
the  following  linear  differential  equation  for  the  joint  density  con- 
ditioned on  the  bearing  0€E[O,2n): 


3p0U,t) 


7 p0(x,t)  - y (x, t)  Pe(x,t) 


The  solution  of  Eq . (4.7)  for  initial  condition  pQ(x,0)  - p„(  ) is 

U 0 


P0(x,t)  - Po(x-v0t)  exp  / y[x- v0(t-s),s]d: 


■ I 

for  x£E‘  and  tf.[0,»).  Tius,  the  probability  of  detect  ing  a fleeing 


datum  bv  time  te [0,»)  for  initial  location  dens  y p^(x)  and  search 
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if: 


density  y(x,t)  Is  given  by 


where  ■ u(cos  6,  sin  0)  is  the  randomly  oriented  velocity  vector  with 
known  magnitude  u.  This  expression  is  used  extensively  for  the  numerical 
optimization  procedure  of  Sec.  5.0. 

4.4  COMPARISON  WITH  A RANDOM  SEARCH  MODEL 

Prior  to  the  work  of  Heilman  [31  and  the  formulation  presented 
here,  problems  concerning  the  search  for  a moving  target  were  often 

♦’rnn#' K \t  ti'in  nf  M/~\  n t-  n C'  o r 1 o clmn!  if  lr»«c  r\  — f Krv  rin/liiin  ondVeh  mo  A ft  1 
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of  Koopman  [ 1 ) . Simulations  Involve  relatively  high  program  develop- 
ment and  execution  costs,  whereas  the  random  search  model  omits  many 
important  features  of  the  operational  search  problem,  such  as  the 
searcher's  actual  path  and  realizable  target  morion.  After  a brief  dis- 
cussion of  the  random  search  model,  this  subsection  solves  a simple 
fleeing— datum  search  problem  using  that  model  and  compares  the  result 
to  the  solution  obtained  using  Eq.  (4.9). 

The  random  search  model  is  predicated  on  the  following  assumptions: 

1.  The  target’s  position  is  uniformly  distributed  in  a region 
of  area  A and  maintains  that  distribution  for  all  te(0,®). 

2.  The  search  path  is  random  in  A in  the  sense  that  disjoint 
sect'ons  of  the  path  are  distributed  uniformly  and  in- 


dependently in  A. 
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3.  The  searcher's  law  of  detection  is  such  that  detection  is 
certain  within  range  d of  the  target  and  cannot  occur 
outside  d. 

Assumptions  1 and  2 do  not  correspond  well  to  physical  imperatives  but 

it  is  assetted  [1]  that,  in  situations  where  target  and  searcher  are 

moving  in  complex  paths  at  varying  speeds,  the  assumptions  arc  reasonable. 

Assumption  3 is  the  well-known  definite-range  law  of  detection,  which 

is  often  used  as  an  approximation  for  a device  with  lateral  range  curve 

00 

q(£).  In  this  case,  d is  taken  to  be  £ q(£)d£. 

Assumptions  1 through  3 lead  to  the  well-known  random  search  formula 
for  the  probability  of  detection 


P - 1 - exp 


where  L is  the  length  of  the  searcher's  path  in  A.  If  we  now  consider  the 
search  area  to  be  a circle  with  Increasing  radius  R(x)  * + ut  (l.e., 

the  radius  is  increasing  at  the  rate  at  which  a fleeing  datum  is  assumed 
to  be  moving  outward),  then  the  probability  of  detection  may  be  approxi- 
mated by  the  random  search  formula  for  an  expanding  area  [6] 

P(r)  - 1 - exp^-  nRJ^d+-T-)-j  • <4-10> 

where  w is  the  searcher's  speed  and  the  search  commences  at  t ■ 0.  We 
will  use  this  expression  as  representative  of  the  random  search  formula- 
tion for  a f leeing- datum  search  problem. 


\ 
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In  order  to  use  Eq.  (4.9),  we  must  first  understand  the  tmpli  -ationa 
of  using  the  definite-range  law  of  detection.  Strictly  speaking,  it  re- 
sults in  an  unbounded  and  discontinuous  search  density  y(x,t),  thereby 
violating  assumptions  (3.11)  and  (3.12).  Although  a formal  argument 
employing  the  limit  of  a sequence  of  bounded  and  continuous  functions 
could  be  made,  a more  direct  way  to  justify  the  use  of  the  definite- 

range  concept  comes  from  a physical  interpretation  of  Eq.  (4.8). 

2 

For  fixed  x££  , 6e[0,2n),  and  tG{0,°°),  Eq.  (4.8)  weighs 
exp[-  I y[x-Vg(t-s),£»] dsj  by  the  value  of  the  initial  density  at  x-v^t. 
For  search  path  Zj(s)£E^,  s€i{0,®),  Y<y,s)  * YjJy.z^s)]  Is  Infinite 
when  |y-Zj(s)|  ^ d,  and  zero  otherwise-  The  motion  of  the  target  im- 
plicit in  the  integral  ^ y1 (x  - (t-s) , is  a straight-line  path 

from  y^  ■ x-v^t  at  s * fl,  to  y.  ■ x at  s * t.  Therefore,  If  the  target 

moving  on  this  path  ever  comes  within  range  d of  the  searcher  at  z^s), 
t 

then  f y.(x-v  (t  s),  z (s)Jds  - <*>,  and  p (x,t)  - 0;  otherwise, 

si  10  1 0 

jT  Yj  (x  - v0  (t-a) , z,(s)]ds  - 0,  and  p0(x,  t)  - PQ(x-v0t).  These  observa- 
tions lead  to  a remarkably  simple  numerical  procedure  for  evaluating 
Eq.  (4.9)  when  a definite-range  law  of  detection  is  assumed. 


Consider  the  following  inputs  to  a hypothetical  submarine  search 


problem: 


1.  The  initial  target  density  is  given  by 


p (x)  « ^—7  exp  [-  (x2  + x2)/2o2]  , x » (x  x )«=K2 

0 2wo2  12  12 


where  0 ■ 10  nmi. 
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7,  The  submarine  moves  as  a fleeing  datum  with  speed  u * 3 kr 

3.  The  search  is  undertaken  at  T.  = 3 hr  by  a single  searcher 
using  a def inite- range  law  sensor  with  d ••  2 nmi . 

4.  The  searcher  flies  at  constant  speed  v * 200  kt  along  the 

path  ^(t)  - , z,2(r)  ] Gt'2,  t - t - Tj  , where 

zu(T>  ” r(r)  cosU  [D-r(r>]| 


z 


12 


(t 


\ 

J 


r (t)  sin 


9 


r (t) 


('  - Vf 


and  D = 32  nmi.  Tne  search  terminates  at  = 7 hr. 

The  search  path  defined  in  item  4 is  an  inward  spiral  starting  at 
D = 32  nmi  from  the  origin  when  t = = 3 hr-  and  ending  at  approxi- 

mately 2 nmi  from  the  origin  when  t - Tp  =■  7 hr. 

The  probability  of  detection  for  this  problem  was  calculated  with 
the  help  of  Eq.  .4.9)  and  plotted  in  Fig.  4.1  under  the  heading  "inward 

spiral  search  path."  The  calculation  was  repeated  lor  a second  search 

* 

path  z^(a)  - /.jdp-r),  t G [U.  T2  _ T^  ) , which  is  simply  an  unfolding  of 
the  first  path.  Figure  4.1  plots  the  result  if  the  second  calculation 
under  the  heading  "outward  spiral  search  path."  The  results  of  these 
tvo  calculations  are  Intuitive-:  The  outwatj  spiral  path  quickly  resu'es 


Fig..  4.1 — Probability  of  detection  for  tvo  search  paths  and  a random  search  aodei 
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.t.n  a relatively  high  probability  of  detection,  since  the  searcher  is  in 
a relatively  "dense"  target  area  fcr  about  the  first  two  hours  of  search; 
the  Inward  spiral  path  takes  about  two  hours  before  rapidly  increasing 
the  probability  of  detection,  since,  after  that  time,  the  searcher 
finally  begins  to  move  through  a relatively  dense  target  area. 

Figure  4.1  also  graphs  the  probability  of  detection  calculated 
using  Eq.  (4.10}  for  random  search  of  an  expanding  area.  The  input 
values  for  this  calculation  are  taken  from  the  first  two  calculations 
as  applicable,  i.e.,  u * 3 kt,  d * 2 nmi,  and  w * 200  kt-  Since  the 
search  commences  at  t » T1  * 3 hr,  we  chose  * 1 . ho  + uT^  = 24  nmi. 

This  choice  is  arbitrary,  but  the  main  conclusion  to  be  drawn  from 
Fig.  4.1  does  not  depend  on  it:  the  path-invariai  t random  search  formula- 

tion is  not  a good  estimate  of  search  system  performance  for  known  search 
tactics. 
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5.0  NUMERICAL  OPTIMIZATION 


The  theory  of  optimal  search  for  a stationary  target  has  received 
considerable  attention;  the  recent  book  by  Stone  [2]  presents  a unified 
treatment  for  much  of  that  work.  Saretsalo  [7],  working  with  Heilman's 
quation  for  detecting  a Markovian  target,  has  presented  a necessary 
condition  for  the  optimality  of  the  search  density  y(x,t)  in  the  sense 
of  maximizing  the  probability  of  detection  P(T)  for  fixed  Tg[0,»). 
Saretsalo's  result  is  directly  applicable  to  the  search  formulation  of 
Sec.  3.0,  but  the  question  of  sufficiency  limits  its  utility. 

The  computational  simplicity  of  Eq.  (4.9)  suggests  the  use  of 
numerical  techniques  to  optimize  the  search  for  a fleeing  datum.  Tliete- 
fore,  this  section  considers  the  numerical  optimization  of  a stationary 
search  (i.e.,  a search  with  density  y(x,t)  « y(  .)  for  ail  x£X  and 
t<=  [T  .Tj))  for  a fleeing  datum  with  a circular  normal  initial  density. 

We  constrain  the  search  density  y(x)  to  be  the  result  of  using  n » 4, 

8,  or  12  definite-range-law  sensors  equally  spaced  on  the  circumference 

2 

of  a circle  with  radius  R centered  at  the  origin  of  E . R is  then  chosen 
to  maximize  P(?2),  the  probability  of  detection  over  the  search  interval 
^l'^2^*  Perf°rmance  the  best  circular  pattern  using  eight  sen- 

sors is  compared  to  tha‘  of  an  optimal  square  n.-it.tern  with  nine  sensors 
and  shown  to  be  superior.  Finally,  we  determ  ie  the  effect  of  using  a 
fleeing  datum  circular  pattern  to  find  a target  that  in  fact  randomly 
chooses  a new  course  every  AT  time  unite.  The  -ectior.  concludes  with 
suggestions  for  broadening  o r computational  experience  and  extending 
the  analytical  development  of  the  search  formulation. 
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5.1  STATIONARY  SEARCH  FOB  A FLEEING  DATUM 

This  subsection  performs  a numerical  optimization  of  a constrained 
stationary  search  for  a fleeing  datum.  Initially,  the  target  is  assumed 
to  he  located  according  to  the  circular  normal  density 

PQ(x)  " exP  t “<x^  +Xj)/202]  , x-(xlfx2)GE2  • (j.1) 

2iro 


Target  motion  for  t _>  0 is  given  by  the  f leeing-datum  transition  prob- 
ability function  of  Eq.  (4.6).  We  restrict  our  attention  to  the  class 
of  search  densities  defined  by 


> .x,t)  - / y(x) 


te  (o,^) 


te.  I‘i1»t2) 


te  (t2,~) 


(5.2) 


A search  using  a member  of  this  class  of  densities  is  referred  to  as  a 
stationary  search  with  time  late  and  duration  T2~T^.  Tlie 
numerical  optimization  seeks  to  maximize  P(T2),  the  probability  of 
detection  during  the  interval  over  a specific  family  of  sensor 

configurations  that  satisfy  Eq.  (5.2). 

The  solution  of  Eq.  (4.7),  the  fleeing- datum  linear  search  equation 
conditioned  on  target  heading  8£[0,2ir),  when  y(x,t)  is  given  by  Eq. 
(5.2),  is 


(x- V) 


k0 

,(X'  * J Pq(x  - v^t)  expj-  J y[x  - Vg (t-s)  ]dsj 

( 


te  [o,t1> 


telTrT2) 
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where  v^  “ u(cos  6,  ein  6)  is  the  target  velocity  vector  with  known 
speed  u.  Hie  probability  of  detection  in  is  thus  given  by 


2n  r- 

P(T2)  " 1~2*  J f PO^X_V0T2)  exp  ' J 

x e-o  L r. 


y{x  - v0(T2~8) ]ds |d9dx 


i 


(5.3) 

where,  for  the  numerical  work  below,  we  assume  that  p^(x)  is  circular 
normal . 

The  search  is  assumed  involve  the  use  of  n defir.itt -range-law 
sensors  with  detection  range  d.  The  sensors  are  assumed  to  be  inde- 
pendent in  the  sense  of  Eq.  (4.1).  (Subsection  4.4  discusses  the  compu- 
tational implications  of  assuming  definite- range-law  sensors.)  We  seek 
the  optimum  placement  of  the  n sensors  such  that  the  probability  of 
detection  calculated  from  Eq.  (5.3)  is  maximized.  However,  we  do  not 
attempt  a global  optimization,  but  instead  constrain  the  sensors  to  lie 
equally  spaced  on  the  circumference  of  a circle  centered  about  x • 0, 
the  mean  of  p^(x) — see  fig-  5.1  for  n * 8 sensors — and  then  calculate 
the  radius  R of  the  circle  that  maximizes  PfT^)  As  long  as  the  cir- 
cular detection  areas  of  the  individual  sensors  do  not  overlap  (as  will 
be  the  case  for  our  analysis),  the  circular  pattern  is  reasonable  in  view 
of  the  radia]  symmetry  of  both  the  initial  density  and  the  fleeing  datum 
motion. 

Placing  the  optimization  in  the  context  of  a search  for  a fleeing 
submarine  by  use  of  moored  acoustic  sensors  (i.e.,  sonobuoysl,  the 
following  values  were  chosen  for  the  input  parameters:  n = 4,  8,  or 

i2  sonobuoys;  u 1 6 or  12  kt;  cj  = 20  or  40  nm.: ; and  d = 4 or  10  nmi. 


■a 


A 


•J 
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The  tine  late  and  search  duration  ••  T were  fixed  at  4 and  6 hr, 
respectively.  The  table  below  lists  the  radius  of  the  optimum  circular 
sonobuoy  pattern,  and  resulting  probabilitv  of  detection,  for  each  com- 
bination of  the  parameters  n,  u,  o,  and  d.  Figure  5.2  plots  p(x,t), 
t ■ 7 hr,  the  conditional  density  midway  into  the  search  for  case  14  in 
the  table.  The  darker  an  area,  the  more  likely  the  target  is  within, 
that  area,  assuming  an  unsuccessful  search.  To  show  the  sensitivity  of 
PfT^) , the  detection  probability,  to  changes  in  R,  the  radius  of  circular 
sonobuoy  pattern.  Fig.  5.3  plots  ^^2)  against  R for  cases  10,  12,  14, 
and  If . 

* 

Figure  5.4  plots  P (T2),  the  probability  of  detection  for  the 
optimum  circular  sonobuoy  pattern,  against  n,  the  number  of  sonobuoys 
deployed.  Each  of  the  eight  '-urves  corresponds  to  a unique  combination 
of  u,  o,  and  d.  The  plot  could  be  used,  say,  to  conclude  that  a search 
for  a fleeing  submarine  using  only  four  sonobuoys,  each  with  a 10  nmi 
detection  range,  is  about  as  effective  as  a search  using  12  sonobuoys, 
e«  ch  having  a 4 nmi  detection  range. 

To  demonstrate  the  effectiveness  of  the  circular  sonobuoy  pattern, 

we  calculated  the  probability  of  detection  for  a square  pattern,  using 

nine  sonobuoys  as  shown  in  Fig,  5.5.  We  chose  the  distance  between 

sonobuoys  L to  maximize  P(T2)  for  speed  u * 12  kt,  detection  range 

d * 10  nmi,  and  initial  standard  deviation  0 s*  20  nmi — as  in  all  cases 

using  the  circular  patterns  • 4 hr  and  T2  *■  10  hr.  The  maximum 
* * 

P (Tq)  - 0.35  was  achieved  for  I,  =»  53  nmi;  Fig.  5.6  plots  p(x,t), 

it 

* 

t « 7 hr  for  L « L . Thus,  the  optimal  square  pattern  with  mne  sono-- 
buoys  was  slightly  less  effective  than  the  optimal  circular  pattern 
using  eight  sonobuoys  (case  8) . 
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(1)  At  t - 0,  target  has  circular 
normal  distribution 
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Fig.  5.3 — Probability  of  detection  versus  radius  for  representative  circular 

sonobuoy  patterns 


Notes:  (1)  At  t * 0,  target  has  circular 

nc-mal  distribution 

(2)  Search  is  from  T.  *•  4 hr 
to  T2  - 10  hr  1 
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^8*  5,4-  Probability  of  detection  for  optimum  circular  sonobuoy 
patterns  used  to  detect  a fleeing  submarine 


Notes:  Target  speed:  u = 12  kt 

Initial  standard  deviation:  a - 20  nmi 

Sonobuoy  detection  <"ange:  d = 10  nir.i 

Distance  between  sonobuoys:  L = 53  rinii 

Time  late:  T^  = 4 hr 

Density  plot:  p(x,  t),  t 7 hr 

!'tg.  5.6  Target  density  p(x,t.)  tor  .m  unsuccessful  search  using  a 


square  pattern  with  n " 9 sonobuoys 
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5,2  EFFECT  OF  TARGET  COURSE  CHANGES  ON  OPTIMIZATION 


We  now  consider  the  effect  of  deploying  an  optimal  circular  sono- 
buoy  pattern  to  detect  a target  that  is  assumed  to  move  as  a fleeing 
datum  when  in  fact  it  randomly  chooses  a new  course  every  AT  time  units. 
Recall  that  a f ceing  datum  moves  on  a straight-line  path  after  a single 
random  heading  election  at  t = 0.  As  in  subsection  5.1,  we  assume  that 
the  target  distribution  at  t = 0 is  circular  normal,  and  that  the  search 
commences  at  time  late  > 0,  The  target  is  assumed  to  move  as  a 
fleeing  datum  until  t = ; however,  at  t = T^,  the  target  selects  a 

new  course  from  a uniform  distribution  on  f0,2n).  Furthermore,  the 
target  continues  to  select  a new  course  every  AT  time  units  corresponding 


fn  ^ ony^  ^ tlhdt  fcCCCTSCG 


uc^iu^iucii  L dl  L 


and  subsequently  attempts  to  avoid  detection. 

At  t = T , the  target  location  density  is  the  solution  of  the 
ssarch-free  equacion  (3.32)  for  a fleeing  datum.  Specifically,  for 
x (x^ , x^) G E ' , 


i 


pU.Tj)  = — ^ exp  [ -;>-2+u2tJ)/2o2  ]Io  (ru'iyc/)  , (5. A) 

2ro 

2 2 2 

where  r - x.  -{-  x„,  and  I denotes  the  ordinary  Bessel  function  of  order 
1 2 o 

zero  with  an  imaginary  argument.  Since  the  search  commences  at  the  tine 
of  the  first  new  course  selection,  Eq.  (5.4)  can  be  interpreted  as  the 
initial  density  of  a f leeing-dat  ur.:  search  problem  with  time  late  zero 
and  search  ducat io  AT.  The  joint  density  at  t = + AT  can  then  be 

computed  using  !iq.  (4.8)  as 
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P(x,T^  + AT) 


2tt 

™Il 


P(x-v  AT,T  ) 
81  1 


^-hvr 

exp[v  1 


y(x  - va  (T  + AT-s' 
®1  1 


j dsj u9 


The  probability  of  no  detection  in  the  interval  JT^,T^+AT)  is 
simply 

QAT(1)  *y*pfx,Ti  f AT)dx 

From  Eq.  (3.28),  the  conditional  density  at  t ” + AT  is 

P (x, Tx  + AT)  - Q~^(l)  p(x,Tj  + AT) 

Now  u9ing  p(x,T^+AT)  as  the  initial  density  for  a fleeing  datum  search 
problem  in  the  interval  [T.  + ATyf^  /AT)  , we  cap.  compute  0^(2) , the 
probability  of  no  detection  in  that  interval,  and  so  on.  This  procedure 
obtains  the  probability  of  detection  in  the  interval  10, T^)  given  by 

m 

w - 1 • n Q«i<k>  • <5-5> 

k”l 


where  AT  is  conveniently  chosen  such  that  mAT  = - T^ . 

For  case  10  in  the  table,  the  radius  of  the  optimum  circular  pattern 
* 

of  eight  sonobuoys  is  R * 32  nmi.  Holding  all  other  inputs  the  same 
as  in  case  10  (i.e.,  submarine  speed  u = 6 kt,  initial  standard  devia- 
tion a * 20  nmi,  and  sonobuoy  detection  range  d ® 10  nmi),  we  use  this 
circular  pattern  to  compute  P (T,,)  from  Eq.  (3.5)  as  a function  of  id,  the 
number  of  course  changes  in  the  6 hr  interval  (T.,,^).  The  results  are 
graphed  in  fig.  5.7.  The  value  of  at  to  = 0 is,  of  course,  the 
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probability  of  detection  for  the  optimum  circular  sonobuoy  pattern  when 
the  submarine  moves  as  a simi  Le  fleeing  datum.  It  is  of  some  interest 
to  note  that  P^d^)  > Pq(T2)  =0.58  for  m * 1 , 2, . . . , 7 . A submarine 
attempting  to  avoid  detection  by  a circular  sonobuoy  pattern  optimized 
tor  a fleeing  datum  would  thus  actually  decrease  its  probability  of 
escape  by  making  a few  additional  course  selections.  Fun  lermore.  Fig. 

5.7  suggests  that  an  optimum  f leeing-da turn  search  pattern  is  not  signi- 
ficantly compromised  by  any  number  of  course  changes  in  the  inter- 
val [T^.Tj)  (i.e.,  |p  (T2>  - P0(I2)  | < 0.03  for  all  m>7).  Although  our 
example  is  somewhat  idealized,  these  results  may  warrant  further  in- 
vestigation with  more  realistic  inputs. 

As  a final  calculation  tor  the  inputs  of  case  10  and  m ~ 6 course 

changes.  Fig.  5.8  plots  P (T,.)  from  Eq.  (5.5)  against  R,  the  radius  of 

ml 

a circular  sonobuoy  pattern.  The  optimum  circular  pattern  of  eight 
sonobuoys  is  shown  to  have  a radius  of  28  nmi  when  the  submarine  is 
assumed  t > select  a new  bearing  every  hour.  This  compares  \ Lth  an 
optimum  radiut  of  32  nmi  for  the  simple  fleeing  datum  motion. 

5 ■ j EXTENSIONS 

The  purpose  of  this  report  is  to  1)  render  tne  search  formulation 
accessible  to  operations  analysts  without  sacrificing  mathematical 
r i;cr;  2.)  narrow  the  gap  between  certain  analytical  assumptions  and  the 
behavior  of  actual  targets  and  existing  search  devices;  and  3)  illustrate 
the  utility  of  the  formulation  with  simple  numerical  op  ■:  Inization  pro- 
cedures. The  developments  herein  should  be  regarded  as  an  introduction  to 
this  work.  Son-.-  suggestions  ft  r further  research  ere  given  below 
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Fig. 
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Notes:  (1)  Search  Inputs: 

•Number  of  sonohuoys:  n » 8 
•Sonobuoy  detection  range:  d ■ 10  nmi 

•Time  late:  Tj  ■ 4 hr 

•Search  duration:  r?  - Tj  » 6 hr 

(2)  Submarine  inputs: 

•Standard  deviation  at  t * 0:  o • 20  nmi 
•Speed:  u * 6 kt 

•First  course  change  at  1 4 hr 
•Later  course  changes  every  hour 
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Radius  of  sonobuoy  pattern,  R (nmi) 


5.8 — Effect  of  sonobuoy  pattern  radius  on  probability  of  detection 
for  a submarine  that  changes  course  every  hour 
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The  numerical  analysis  of  this  section  would  be  taore  relevant  to  an 
actual  search  If  a modest  effort  were  made  to  Improve  the  characterization 
of  target  and  searcher  behavior.  For  example,  by  obtaining  analytical 
solutions  of  F,q.  (3.30)  for  1)  the  class  of  conditionally  deterministic 
tar  ;et  motion  or  2)  more  general  diffusion  motions,  we  could  pet  form  the 
numerical  work  of  this  section  for  a much  broader  and  realistic  set  of 
targets . 

With  regard  to  the  search  process,  it  would  be  desirable  to  solve 
numerical  problems  foj  search  devices  with  more  sophisticated  laws  of 
detection  than  the  definite-range  concept.  Perhaps  an  integration -type 
sensor  (see  subsection  4.2)  subject  to  "convergence  zone"  phenomena 
| tl J would  be  a more  realistic  device  to  study.  More  generally,  future 
numerical  work  should  include  the  constrained  optimization  of  a moving 
searcher  problem — i.e.,  from  a given  family  of  paths  that  a moving 
searcher  might  follow,  find  the  path  that  results  in  the  highest  prob- 
ability of  detection  In  a given  search  problem. 

On  a more  fundamental  level,  work  should  be  done  on  the  problems 
of  evasive  targets,  false  targets  and  false  contacts,  multiple  targets, 
and  search  devices  chat  detect  as  a function  of  the  relative  orientation 
of  target  and  searcher.  Finally,  for  a general  Markoviar.  search,  we 
have  yet  to  solve  the  problems  of  obtaining  sufficient  conditions  for 
an  optimal  search  density  and  hew  to  obtain  such  a density  having  found 
those  conditions. 
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